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ABSTRACT 

We examine the relation between oxygen abundances in the narrow-line regions 
(NLRs) of active galactic nuclei (AGNs) estimated from t he optical emission lines 
throu gh the strong-line method (the theoretical calibration of IStorchi-Bergmann et all 
[T^ , via the direct Tg-method, and the central intersect abundances in the host galax¬ 
ies determined from the radial abundance gradients. We found that the Tg-method 
underestimates the oxygen abundances by up to ^2 dex (with average value of ~ 0.8 
dex) compared to the abundances derived through the strong-line method. This con¬ 
firms the existence of the so-called “temperature problem” in AGNs. We also found 
that the abundances in the centres of galaxies obtained from their spectra trough 
the strong-line method are close to or slightly lower than the central intersect abun¬ 
dances estimated from the radial abundance gradient both in AGNs and Star-forming 
galaxies. The oxygen abundance of the NLR is usually lower than the maximum at¬ 
tainable abundance in galaxies (^2 times the solar value). This suggests that there is 
no extraordinary chemical enrichment of the NLRs of AGNs. 

Key words: galaxies: general - galaxies: evolution - galaxies: abundances - galaxies: 
formation- galaxies: ISM 


1 INTRODUCTION 


Active galactic nuclei (AGNs) and star-forming regions show 
prominent emission lines of heavy elements that can be eas¬ 
ily measured, even for objects at large redshifts. The inten¬ 
sity of these emission lines depends on the metallicity, which 
makes them convenient tracers of the chemical evolution of 
the Universe. 

The abundance of a given element in AGNs and in 
star-forming regions can be derived from measurements 
of the relative strengths of the emission lines of its ions 
and of the electron temper ature and density of the gas 
llOsterbrock &: Ferlandl [2OO6II . Oxygen has been generally 
used as tracer of the metal l icity (Z) of star-forming regions 

(e-g-, 


l2006a l and of AGNs ( Izotov et al.l 2010l : iGroves et al 


Tremonti et al.l [20041 : Thuan et al.l 20^; Naeaoet^ 


I2OO6I. among others). Methods for the abundance deter¬ 
minations in star-f orming regions have been discussed in 
many papers (e.g., IStasiiisk^boOdh . The consensus is that 
bona-fide determi nations of Z are o nl y obtained by usin g 
the Te-method (iHaeele et al.l I2OO6I : iHagele et al.l I2OO8I I. 
This method is based on the determination of the electron 
temperature (Te) from the ratios between intensities of two 
emission lines originating in transitions from two levels with 
different excitation energies of a same ion, such as the ratios 
[O III] (A4959-bA5007)/A4363, [S iii] 1A9069+A9532)/A6312 
and [Nii](A6548-l-A6584)/A5755 (see lOsterbrock fc Ferlarid 
120061 .1 Unfortunately, these line ratios cannot easily be 
measured in di stant objects and/o r in objects with low 
excitation fe.g.. iBresolin et al.ll2005ll . preventing the use of 
the Te-method. 


E-mail:olidors@uni vap.br 


Therefore, the oxygen abundances in Hll regions are 
us ually estimated us ing t he strong-line meth od pioneered 
by IPagel et al.l lll979h and lAlloin et al.l (Il979ll . The princi- 
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pal idea of the strong-line method is to establish the re¬ 
lation between the (oxygen) abundance in an Hll region 
and some combination of the intensities of strong emis¬ 
sion lines in its spectrum, i.e., the position of the objects 
in the strong-line diagram is calibrated in terms of abun¬ 
dance. Therefore such a relation is often called a “calibra¬ 
tion”. The calibrations of strong emission lines in terms 
of the oxygen abundances can be accomplished by three 
methods: (i) theoreti cal calibrations based o n photoion¬ 
ization models (e.g., Kewlev fc Dopit^ I2OO2I : iDors etliD 
l201ll : [Dors k, CopettabOOSli . (ii) empirical calibrations based 
on H II re gions with ab u ndances derived through the 
Te-method Ipilvu einl 20 0d. 200ll : IPilvuein fc Thuan l2005l : 
IPerez-Montero fc Continil 20091 . among others), and (iii) hy¬ 
brid calibrations where both photoionization models and H ll 
regions with Te-based abundances are used as calibration 
data points JPettini fc Pagell l2004ll . The metallicity scales 
produced by these d ifferent methods may d iffer from each 
other by up 0.7 dex (iKewlev fc Ellisonll2008h . 

Here we use the metallicity scale in Hll regions de¬ 
fined by the abundances determined through the classical 
Te-method. This is because in the Te method the physical 
conditions in the nebulae, which are essential ingredients 
in order to calculate the abundance, are derived directly 
from observations. Otherwise, if the physical conditions and 
consequently the abundances are derived by models they 
will be less certain. This is because of the parameters of 
the model such as the ionizing sourc e, geometry, etc . , tend 
not to be sufficiently realistic (see iDors fc Copettil l2005l : 
iKennicutt et al.l[2003l : IVieeas l[2002l '). 

AGN metallicities are usually e stimated through strong¬ 
line theoretical calibrations (e.g., IStorchi-Bergmann et al.l 
1 19981 ). Even solar or s u perso lar abundances have been 
found. IPichardson et ahl (120141 ) established a sequence of 
photoionization models to reproduce the narrow-line re¬ 
gions (NLRs) of AGN spectra and found that models 
with a metallicity of 1.4 Zq prov ide the best agree - 
ment with the observational data. iGroves et al.l (l2006l ) 
used the photoionization models to analyse the emission 
lines in the spectra of the NLRs from the Sloan Digi¬ 
tal Sky Survey (see iTremonti et al.l l2004l) and found su- 
persolar metallicities for typical Seyfer t galaxies (see also 
Batra^_&_^aldwm[_[20l4jDu/et_alJj201^JWMig_^£^jJ_l201l|; 


Dhanda^_et_al _ 200^^ Baldwiriet_alJ 200^;_ Hainumie£_aJj 


200^;_ StorchTBerCTiannet_a^ 1 19981 : Ferland et al. 19961 : 
Hamann fc Ferlandl Il993l . Il992l ). High metallicities at the 


centres of spiral galaxies are also obtained by the extrap- 
olatio n of radial abundance gradients to the central re¬ 
gions ifVila-Gostas fc Edmunds lll992l:IZaritsky et al. Ill994l: 


Zee et al.lll998l : IPilvugin et al.l 2004 . 2007 : Gusev et ahl 
201 I among others). 


In contrast, low metallicities in AGNs h ave been ob¬ 
taine d when the Te method was used (e.g., IZhang et ahl 
I2OI3I ). The disagreement between the abundances obtained 
through the strong-line methods and through the Te-method 
is the so-called “tempe rature problem”. On the other hand, 
lizotov fc ThuanI (l2008l ) determined the abundances through 
the Te-method in the AGNs in four dwarf galaxies and found 
metallicities between 0.05 and 0.2 of t he solar value , whic h 
is a typical abundance of those galaxies. lAlloin et al. liHii), 
also using the Te-method, estimated the abundance in the 
Seyfert 2 nucleus of the galaxy ESO 138 G1 classified as E/SO 


dLauberdl 19821) and found Z « 0.4 Zq. Thus, the Te method 
seems to produce realistic abundance estimations in the low 
metallicity range, i.e., the temperature problem seems to 
disappear at low metallicities. Finding the correct method 
for the determination of AGN metallicities and establishing 
a relation between the AGN metallicity and the metallicity 
of the host galaxy are very important challenges. It should 
be noted that the AGN metallicity is sometimes adopted 
as a surrogate met allicity of the host galaxy, mainly at 
high redshifts (e.g. , iDors et al.ll2014l : iMatsuoka et al.ll2009l : 
iNagao et al.ll200^ . 

Nowadays, a large number of measured intensities of 
emission lines sensitive to the electron temperature are avail¬ 
able for AGNs in the literature. This also holds for de¬ 
terminations of oxygen gradients in a large sample of spi¬ 
ral galaxies. This provides the possibility for an analy¬ 
sis of the metallicity determinations in AGNs. The main 
goals of the current study are: (i) To investigate the dis¬ 
crepancy between AGN abundances derived through the 
strong-line methods and through the Te method. The 
magnitude of thi s discrepancy i s well known for star¬ 
forming regions (iPilyugi nl l2003l : I Kewlev fc EllisonI l2008l : 
iLopez-Sanchez fc EstebanI 120101 7 among others), however, 
this problem has received less attention in the case of the 
NLRs of AGNs. (ii) To compare the abundances of the NLRs 
of AGNs and nuclear star-forming regions estimated from 
their spectra and central oxygen abundances of the host 
galaxies derived from the radial abundance gradients. Dif¬ 
ferences between these values could be evidence in favour 
of the existence of gas infall onto the center of a galaxy or 
extraordinary chemical enrichment of these objects. 

The present study is organised as follows. In Section 
2 we compare the NLR abundances determined through 
the Te-method and the stron g emission-line relations of 
IStorchi-Bergmann et ah I ll 19981. hereafter SB98). We exam¬ 
ine the relation between the abundances in the NLRs esti¬ 
mated from their optical emission-line spectra and the cen¬ 
tral intersect abundances in the host galaxies determined 
from the radial abundance gradients in Section 3. The con¬ 
clusions are given in Section 4. 


2 ABUNDANCE DETERMINATIONS IN 
NLRS: Tb-METHOD VS. RELATIONS BY 
SB98 

2.1 AGN sample 

Intensities of narrow emission lines of AGNs were compiled 
from the literature. Our selection criterion was the presence 
of flux measurements of the narrow optical emission lines 
[Oii]A3727, [Oiii]A4363, H,5, [Oiii]A5007, Ha, [Nii]A6584, 
and [S ii] AA6717,6731. We only considered galaxies whose 
nuclei were classified as Seyfert 2 (Sy2) and Syl.9. Sequences 
of photoionization models, not including the shock gas but 
only a power law radiation of the ionizing source, are able 
to reproduce measurements of strong narrow e mission lines 
of Sy2 and Syl.9 for a large spectral range (e.g. jRiffel et ahl 
l2013al : [Pors et ^l2012l : iGroves et akllioo^ ). Hence, the se¬ 
lection of these AGN types minimises the effects of shock 
gas excitations and ionizations, which are not considered in 
our abundance determinations. 
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We list in Table [T] the identification of the objects, the 
emission-line intensities (relative to H/3=1.0), the spectral 
classification taken from the neeQ, and bibliographic refer¬ 
ences from which the data were taken. The sample consists 
of 47 local AGNs (redshift 2 < 0.1) observed with long- 
slit spectroscopy. In the cases where the reddening correc¬ 
tion was not performed in the original works (indicated in 
Table. [T|), we dereddened the emission-lines comparing the 
observed Hq/H/ 3 ratio w ith the theoretical value of 2.86 
dHummer fc Storevlll997tl . obtained for an electron temper¬ 
ature of 10 000 K and an electron density of 100 cm“®. 

The emission-line flux errors are reported in the orig¬ 
inal works for only 4 out of the 47 objects of our sample. 
Therefore the errors are not given in Table [T] Typical er¬ 
rors of the emission-line intensi ties are about 10 to 20 per 
cent (e.g., iKraemer et al. |[T99i), which yield uncertainties 
in the ox ygen abundance estimations of abou t 0.1 to 0.2 


dex ( e.g., Rosa^_eLalJ 


I 2 OIII : iKennicutt et al.l 


201T; iHagele et al.l I2OO8I : iDors et H] 


20031) . Hereafter we will assume that 


the oxygen estimations from the Te-method and the rela¬ 
tions by SB98 have an uncertainty of 0.15 dex, an average 
of the values above. For the objects which the flux of the 
line [Olll]A4959 is not available, this was calculated from 
the theoretical relation /[O lll]A4959=/[0 lll]A5007/3.0. 


2.2 Relations of SB98 

SB98 carried out NLR mod el calculations usin g the pho¬ 
toionization code Cloudy (iFerland et al.l Il996ll and sug¬ 
gested two relations for the abundance determinations in 
the NLRs of AGNs. The first one is 

= 8.34-f (0.212*)-(0.012 *2)-(0.021/) 

-k (0.007 xy) - (0.002 x^y) + (6.52 x lO"'^ y^) 
+ (2.27 X lO-'^ xy^) + (8.87 x 10“® x^y^), 

where * = [N ii]AA6548,6584/Ha and y = 

[O iii]AA4959,5007/H/l. The second relation of SB98 
is 

(0/H)gggg 2 = 8.643 - 0.275 u-k 0.164^2 

-k 0.655 u - 0.154 uu - 0.021 (2) 

-k 0.288u^ -k 0.162 mu^ -k 0.0353u2w^ 

where u = log([0 ii]AA3727,3729/[0 iii]AA4959,5007) and v 
= log([N ii]AA6548,6584/Hq). Both calibrations are valid for 
8.4 ^ 12-klog(0/H) S 9.4. 

The dependence of these relations on the density should 
be taken into account. This dependence is given by the ex¬ 
pression considered by SB98 

(0/H)fi„,i = 0/H - 0.1 log(Ne/300), (3) 

where is the electron density in cm“® and the correction 
is valid for 100 cm“® < 10“^ cm“^. It should be noted 

that the value of this correction exceeds 0.1 dex for high- 
density objects only, i.e., for objects with electron densities 
Ne >3x10® cm"®. 


^ The NASA/IPAC Extragalactic Database (ned) is operated by 
the Jet Propulsion Laboratory, California Institute of Technology, 
under contract with the National Aeronautics and Space Admin¬ 
istration. 


2.3 Te-method 


We use the emission-line intensities listed in Table [T] 
to determine the oxygen abundances (0 /H)t^ of the 
narrow-line regions through the classic Te-method. We wil l 
follow the methodology describ ed in iDors et al.l (l2013f) 
based o n the equatio n s give n by |Percz:Montero_fc_Xkgi^^ 
ll2009ll. Hagele et al. (l200^'l. IPerez-Montero et al.l ll2007h . 
and Perez-Montero fc DfazI ll2003h . 

The electron temperature in the high ionization zone 
(referred to as fa) for each object was obtained from the ob¬ 
served line-intensity ratio i?03=[0 iii](A4959+A5007)/A4363 
using the expression 


t3 = 0.8254 - 0.0002415Ro3 + 


47.77 

Ro3 


(4) 


with ts in units of lO'^K. The electron temperature of the 
low ionization zone (referred to as 12 ) was derived from the 
theoretical relation: 


= 0.693 -k 0.281. 


(5) 


The 0» and 0> ionic abundances were computed through 
the relations: 

,0 , ,0++, , ,7(4959)-k 7(5007). 

12 + log( —) = log(^-+6.144 


lim 

_^L2W _ 0.55 log fa. 


and 


7a 

, ,0+, , ,7(3727) 

12 + log(jj^) - log( 


(6) 


) + 5.992 


1 583 

+ -0.681 logt2 + log(l + 2.3ne).(7) 

C2 

(1) Finally, we assumed 

O _ 0 + 0 ++ 

H - IT + 

for the determination of the total abundance. 


2.4 Abundance results 

Electron densities (Ae), Olll electron temperatures (Te), 
and the oxygen abundances estimated using three differ¬ 
ent ways for the objects in our sample are reported in Ta¬ 
ble O The oxygen abundances of some objects can not be 
determined through the relations given by SB98 because for 
those objects the estimated abundances or electron densi¬ 
ties are outside of the validity ranges considered by these 
authors, i.e., beyond of the range of values for which the 
relations are defined. The electron densities, Ne, for the ob¬ 
jects in our sample were computed from the line intensity 
ratio [S ii]A6716/A6731 using the temden routine of the neb¬ 
ular package of IRAI0. Most of our objects have values in the 
range 0.01 < rie < 0.12, where rie = Ae/(10'^cm“®). These 
values are somewhat higher than the ones der ived for H ll re¬ 
gions , which were found to be rie < 0.06 (see lKrabbe et alTI 
|2014 IConetti et al. l[2000f) . However, they are sufficiently 
low so that one does not need to consider the contribution 


^ Image Reduction and Analysis Facility, distributed by NOAO, 
operated by AURA, Inc., under agreement with the NSF. 
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Table 1. Emission-line intensities of AGNs relative to H/3=1.00 compiled from the literature. 


Object 

[O ii]A3727 

[Oiii]A4363 

[Oiii]A4959 [Oiii]A5007 

Ha 

[Nii]A6584 

[S ii]A6717 

[S ii]A6731 

Type 

Ref. 

NGC 7674 

1.29 

0.12 

3.94 

12.55 

3.70 

3.68 

0.54 

0.64 

Sy2 

1 

IZw92 

2.63 

0.32 

3.52 

10.12 

2.42 

0.97 

0.37 

0.40 

Sy2 

1 

NGC 3393(“b(6) 

2.41 

0.14 

5.47 

16.42 

2.78 

4.50 

0.64 

0.89 

Sy2 

2 

Mrk 176 

3.54 

0.32 

4.69 

14.36 

2.81 

2.99 

0.56 

0.54 

Sy2 

3 

3c33 

4.93 

0.32 

4.10 

12.68 

2.63 

1.76 

0.87 

0.73 

Sy2 

3 

Mrk 3 

3.52 

0.24 

3.99 

12.67 

3.10 

3.18 

0.73 

0.82 

Sy2 

3 

NGC 1068 

1.23 

0.22 

4.11 

12.42 

2.57 

4.55 

0.26 

0.55 

Sy2 

3 

Mrk 573 

2.92 

0.18 

3.89 

12.12 

2.95 

2.47 

0.75 

0.80 

Sy2 

3 

Mrk 78 

4.96 

0.14 

3.88 

11.94 

2.46 

2.32 

0.68 

0.61 

Sy2 

3 

Mrk 34 

3.43 

0.15 

3.68 

11.46 

2.99 

2.18 

0.82 

0.80 

Sy2 

3 

Mrk 1 

2.78 

0.21 

3.56 

10.95 

2.66 

2.21 

0.49 

0.52 

Sy2 

3 

3c433 

6.17 

0.41 

3.27 

9.44 

3.38 

5.13 

1.58 

1.13 

Sy2 

3 

Mrk 270 

5.64 

0.28 

2.96 

8.71 

3.14 

2.93 

1.21 

1.39 

Sy2 

3 

3c452 

4.81 

0.18 

2.40 

6.85 

2.98 

3.58 

1.10 

0.77 

Sy2 

3 

Mrk 198 

2.51 

0.12 

1.79 

5.56 

3.02 

2.26 

0.89 

0.68 

Sy2 

3 

Mrk 268 

3.75 

0.25 

1.55 

4.82 

3.38 

4.94 

1.28 

1.08 

Sy2 

3 

Mrk 273 

8.27 

0.22 

1.52 

4.44 

2.94 

2.62 

0.89 

0.54 

Sy2 

3 

NGC 3227('’) 

3.22 

0.50 

3.57 

10.73 

2.86 

5.01 

1.24 

1.26 

Sy2 

4 

Mrk 6 

2.45 

0.28 

3.37 

10.13 

2.79 

1.79 

0.62 

0.63 

Sy2 

4 

ES0138G1 

2.35 

0.34 

2.93 

8.71 

3.01 

0.68 

0.47 

0.48 

Sy2 

5 

NGC 5643(“)>(f<) 

5.55 

0.54 

4.20 

12.61 

2.66 

2.90 

0.91 

0.71 

Sy2 

6 

NGC 1667 

11.50 

0.59 

2.78 

9.20 

2.80 

6.96 

1.37 

1.17 

Sy2 

7 

Mrk 423 

8.00 

0.40 

2.80 

6.20 

2.90 

3.60 

1.20 

1.00 

Syl.9 

8 

Mrk 609 

1.80 

0.50 

1.80 

5.00 

2.80 

2.60 

0.60 

0.50 

Syl.9 

8 

Mrk226SW 

5.23 

0.08 

1.50 

4.50 

3.35 

2.07 

0.29 

0.25 

Sy2 

9 

NGC 308l(“) 

2.16 

0.23 

4.53 

12.62 

2.73 

2.33 

0.60 

0.62 

Sy2 

10 

NGC 328l(“) 

2.33 

0.22 

2.77 

7.59 

2.64 

2.54 

0.53 

0.60 

Sy2 

10 

NGC 3982(“) 

3.92 

0.03 

5.64 

18.68 

2.76 

2.57 

0.76 

0.82 

Sy2 

10 

NGC 4388(“) 

2.68 

0.15 

3.67 

10.63 

2.71 

1.44 

0.68 

0.60 

Sy2 

10 

NGC5135(“) 

2.01 

0.10 

1.40 

4.47 

2.64 

2.35 

0.37 

0.35 

Sy2 

10 

NGC 5643(“) 

5.11 

0.42 

4.56 

15.4 

2.64 

3.07 

0.97 

0.90 

Sy2 

10 

NGC 5728(“) 

3.41 

0.44 

3.70 

10.98 

2.65 

3.71 

0.41 

0.41 

Sy2 

10 

NGC 6300(“) 

15.48 

1.39 

8.32 

23.32 

2.59 

6.62 

1.42 

1.26 

Sy2 

10 

NGC 6890(“) 

2.77 

0.72 

7.48 

20.05 

2.75 

4.26 

0.65 

0.52 

Sy2 

10 

IC 5063(“) 

5.06 

0.28 

3.37 

10.31 

2.77 

2.67 

0.69 

0.60 

Sy2 

10 

IC5135(“) 

4.05 

0.25 

2.07 

6.88 

2.65 

3.30 

0.49 

0.46 

Sy2 

10 

Mrk 744 

2.38 

0.33 

3.18 

8.84 

2.47 

3.62 

2.83 

2.83 

Sy2 

11 

Mrk 1066 

3.34 

0.08 

1.22 

3.84 

2.76 

2.42 

0.51 

0.55 

Sy2 

11 

NGC 5506 

2.84 

0.14 

2.46 

7.69 

2.84 

2.53 

0.92 

0.99 

Syl.9 

12 

NGC 2110 

4.38 

0.26 

1.61 

4.76 

2.66 

3.76 

1.52 

1.38 

Sy2 

12 

NGC 328l(“b{*>) 

4.30 

0.42 

2.44 

7.34 

2.60 

2.60 

0.78 

0.73 

Sy2 

13 

Akn347(“) 

2.98 

0.42 

4.95 

15.01 

2.65 

3.23 

0.75 

0.75 

Sy2 

14 

UM16(“) 

2.90 

0.22 

4.62 

14.00 

2.72 

1.70 

0.45 

0.45 

Sy2 

14 

Mrk533(“) 

1.59 

0.13 

4.03 

12.23 

2.72 

2.72 

0.39 

0.45 

Sy2 

14 

IZw92(“) 

2.60 

0.34 

3.54 

10.14 

2.78 

1.00 

0.40 

0.43 

Sy2 

14 

Mrk612(“) 

1.88 

0.17 

2.99 

9.37 

2.67 

3.60 

0.74 

0.55 

Sy2 

14 

Mrk622(“) 

10.06 

0.03 

1.90 

5.44 

2.47 

2.33 

0.19 

0.14 

Sy2 

14 
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from the collisional de-excitation (lRubinlll989li or the direct 
recombination of the forbidden lines used in the standard 
abundance determinations. 

We found the electron temperature ts to be in the 
range of 6 000 to 50 000 K, with an average value of ~20 000 
K. This value is somewhat higher than the one found by 


IZhang et al.l (l2013l l for Sy2s, who used data from SDSS DR7 
and found an average value of 14 000 K. 

Fig. [T] shows the comparison between oxygen abun¬ 
dances in NLRs determined through the two methods pre¬ 
sented above. Inspection of Fig.[T]reveals that the Te-method 
provides oxygen abundances lower than the ones obtained 
via the relations of SB98 by up to ~2 dex, with an average 
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T able 2. Object name, electron d ensity, electron temperature, and oxygen abundance obtained via the Tg-method and via the calibrations 
of IStorchi-Bergmann et alJ lll998l i (Equations [T] and [2]l for the objects in our sample. Oxygen abundance value estimations outside the 
validity interval of the calibrations were not considered. 


Object 

Ng (cm-3) 

Tg (lO'^K) 

12-blog(0/H) 
Te-method SB98,1 

SB98,2 

NGC 7674 

1148.0 

1.14 

8.49 

8.72 

9.19 

IZw 92 

822.0 

1.95 

7.90 

8.44 

8.67 

NGC 3393 

2022.0 

1.09 

8.69 

— 

9.24 

Mrk 176 

535.0 

1.61 

8.21 

8.84 

9.02 

3c33 

252.0 

1.72 

8.14 

8.64 

8.80 

Mrk 3 

948.0 

1.49 

8.26 

8.75 

8.96 

NGC 1068 

26993.0 

1.44 

8.40 

— 

— 

Mrk 573 

781.0 

1.34 

8.33 

8.65 

8.92 

Mrk 78 

370.0 

1.22 

8.49 

8.74 

8.88 

Mrk 34 

546.0 

1.27 

8.39 

8.61 

8.85 

Mrk 1 

767.0 

1.50 

8.17 

8.63 

8.90 

3c433 

50.0 

2.38 

7.84 

8.96 

9.04 

Mrk 270 

1027.0 

1.97 

7.99 

8.62 

8.75 

3c452 

50.0 

1.76 

7.95 

8.80 

8.95 

Mrk 198 

111.0 

1.59 

7.88 

8.61 

8.84 

Mrk 268 

260.0 

2.68 

7.52 

8.76 

8.93 

Mrk 273 

50.0 

2.60 

7.73 

8.67 

8.72 

NGC 3227 

647.0 

2.49 

7.75 

8.98 

9.15 

Mrk 6 

647.0 

1.81 

7.95 

8.55 

8.84 

ESO 138 G1 

685.0 

2.22 

7.73 

8.40 

8.53 

NGC 5643 

141.0 

2.36 

7.90 

8.86 

8.96 

NGC 1667 

281.0 

3.12 

7.84 

9.15 

9.06 

Mrk 423 

239.0 

3.14 

7.67 

8.73 

8.80 

Mrk 609 

239.0 

4.41 

7.11 

8.62 

8.92 

Mrk226SW 

296.0 

1.45 

8.08 

8.52 

8.60 

NGC 3081 

693.0 

1.84 

7.85 

8.60 

8.91 

NGC 3281 

974.0 

1.86 

7.84 

8.61 

8.90 

NGC 3982 

819.0 

0.68 

9.54 

8.87 

8.99 

NGC 4388 

343 

1.31 

8.30 

8.54 

8.80 

NGC 5135 

492.0 

1.61 

7.78 

8.57 

8.83 

NGC 5643 

451.0 

1.79 

8.17 

8.93 

8.99 

NGC 5728 

606.0 

2.26 

7.84 

8.86 

9.05 

NGC 6300 

360.0 

2.96 

8.11 

— 

9.17 

NGC 6890 

176.0 

2.11 

8.06 

— 

9.35 

IC 5063 

311.0 

1.79 

8.06 

8.51 

8.65 

IC5135 

471.0 

2.12 

7.78 

8.72 

8.90 

Mrk 744 

606.0 

2.16 

7.76 

8.83 

9.10 

Mrk 1066 

839.0 

1.56 

7.90 

8.54 

8.69 

NGC 5506 

809.0 

1.46 

8.09 

8.60 

8.85 

NGC 2110 

395.0 

2.78 

7.54 

8.73 

8.87 

NGC 3281 

471.0 

2.87 

7.61 

8.66 

8.83 

Akn 347 

606.0 

1.82 

8.09 

8.93 

9.11 

UM16 

606.0 

1.37 

8.35 

8.59 

8.87 

Mrk 533 

1046.0 

1.18 

8.44 

8.72 

9.12 

IZw 92 

805.0 

2.02 

7.87 

8.42 

8.65 

Mrk 612 

75.0 

1.46 

8.10 

8.89 

9.24 

Mrk 622 

64.0 

0.97 

8.76 

8.69 

8.73 


value of ~ 0.8 dex, which confirms the existence of the so- 
called “temperature problem” in AGNs. This discrepancy 
is muc h larger than the one fo und for H ll regions. For ex¬ 
ample, iKennicutt et al.l (l2003f i. using spectroscopy data of 
Hll regions in spiral galaxies, found that O/H abundances 
computed from re l ations based on photoionization models of 
iKewlev fc Dopi^ (l2002l l are higher by 0.2 to 0.5 dex than 
the ones from the Tg-method. 

The abundance discrepancy in AGNs could be caused 
by the presence of a secondary heating (ionizing) source 


in addition to the radiation from the gas accretion onto 
the AGN. This secondary source i s probably related to the 
shock. Indeed, IZhang et al.l (l2013ll showed that the strong 
[Olll]A4363 flux in AGNs (and consequently high electron 
temperature) suggests the presence of some supplementary 
energy source(s), w hich could be due to the presence o f 
shock waves (see also lGontini et aDl2012l : rPrieto et al.ll200^ 1. 
More over, the presen ce of fluctuations of electron tempera¬ 
ture dPeimberd 1 1967^ in the gas phase of AGNs could also 
contribute for the discrepancies found in Fig. [T] 
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12 + log(0/H)y 


Figure 1. The bottom panel shows the comparison between the 
oxygen abundance obtained through the Te method and through 
the two calibrations of SB98 as indicated in each plot. The upper 
panel shows the differences between those O/H estimations. 


In order to examine if Sy2 galaxies have a secondary 
heating sonrce, we performed a simple test. In Fig. [21 the in¬ 
tensity of [O iii]A4363/H/3 versus the R 23 parameter defined 
as R 23 = [/([O II]A3727) -h /([O III]A4959 + A5007)]/7(H/3) 
for the objects in our sample and the ones predicted by a 
grid of photoionization models are shown. 

We consider two models to reproduce the observational 
data: (i) the AGN model to describe the NRL spectra, and 
(ii) the SF-I-AGN model to describe composite (star-forming 
regions -|- A GN) spectra. We us e the version 08.00 of the 
Gloudy code (iFerland et al.l ll 996l) to construct NLR models 
similar to the ones used by Dors_et_^ I ll2012l) , but con sid- 
ering the Table-AGN model ( Mathews fc FerlandlllO^ ) as 
the ionizing source. We considered a large range of values 
for the physical parameters: the number of ionizing photons, 
51 ^ logQ(H) ^ 55; electron densities, 1.0 ^ log(Afe) ^ 4.0; 
and metallicities, 0.05 ^ (Z/Zq) ^ 2. The spectrum of the 
SF-I-AGN model is the sum of the predictions of the SF 
model and the AGN m odel. The SF model is constructed 
using the Cloudy code (iFerland et al.l fl9^ assuming the 
ionizing source to be a stellar cluster with an age of 2.5 
Myr and logQ(II) = 52.85 wh ose spectrum is com puted 
with the STARBURST99 code (iLeitherer et al.ll 19991) . The 
metallicity of the gas phase was considered to be solar with 
an electron density Ne — 100 cm“^. These parameters are 
similar to the ones derived in circumnuclear star-formation 
regions observed in two galaxies ( NGC 1097 and NG C 6951) 
containing a Seyfert 2 nucleus bv iDors et al.l (I2OO8I ). 

We can see in Fig.[2]that the ratio [O lll]A4363/H/3 is un¬ 
derpredicted by the AGN and SF-I-AGN models when com¬ 
pared the observations. This line ratio depends strongly on 
the electron temperature, therefore the electron tempera¬ 
tures in AGNs predicted by our models are lower than the 
ones estimated using the Te-method. Thus, this discrepancy 



109(^23) 

Figure 2. The [O III] A4363/H/3 line ratio as a function of the 
R 23 = [7([0 II]A3727) + 7([0 III]A4959 + A5007)]//(H/3) parame¬ 
ter. The left panel shows the predictions of the AGN model. The 
predictions of SF-I-AGN model are plotted in right panel. The 
points are the objects of our sample (see Table (TJ. The typical 
error (error bars not shown) of the emission-line ratios is about 
10 per cent. The solid lines connect the predictions of the models 
with the same metallicity Z, while dashed lines connect models 
with the same ionizing photon number (Q(H)). 


could be attributed to the presence of gas shock waves prop¬ 
agating at supersonic velocities through the NLRs, which 
does to enhance the intensity of [Olll]A4363 producing 
larger electron temperature values and, consequently, low 
(unrealistic) O/H values when using the Te-method (see 
iNagao et al.ll200ll and references therein). 

Due to the result above, in our paper the Te-method is 
no longer used to compute the metallicity of AGNs. 

SB98 pointed out that an averaged value of Z obtained 
from both relations (Eq. [T] and [2| should be used. Instead 
of doing that we use our models to choose which of the re¬ 
lations by SB98 yields more reliable abundances. For that, 
we use the standard O/H - R23 empirica l diagr am. The R23 
parameter was suggested bv jPagel et al.l (1 19791) as indicator 
of the oxygen abundance when the Te-method cannot be ap¬ 
plied. Although R23 has a strong dependence on metallicity, 
it also depends on the ionizati on degree of the nebu lae, which 
should be taken into account (|Pilvuginll200ll . [2000l ). We com¬ 
pare the predictions of our AGN models with the observed 
O/H - R23 diagram for the objects in our sample, where 
the oxygen abundances are estimated using both relations 
of SB98 (see Fig. [3]). Since the AGN and SF+AGN mod¬ 
els predict similar values of R23 in the zone of Fig. [2] where 
our data are located, only the former models are shown in 
Fig. (21 In this figure we can see that a better agreement is 
obtained when the abundances are determined through the 
first relation of SB98. Hereafter this relation will be used. 
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109(^23) 

Figure 3. 12+log(0/H) vs. log(i? 23 )- The points represent oxy¬ 
gen abundances in our objects estimated through the two rela¬ 
tions (see Sect. [Oil proposed by SB98 as indicated in each plot. 
Lines show our photoionization model results for different values 
of electron densities, as indicated, and the logarithm of the ion¬ 
izing photon numbers (log Q(H)=51-54, with a step of 1 dex). 
The arrow indicates the direction in which the logarithm of the 
ionizing photon numbers increases in each model set. 



log( [NII]>u 6584/Ha) 

Figure 4. BPT [O iii]A5007/H/3 vs. [N ii]A6584/Hck diagnostic di¬ 
agram. Filled squares and open circles stand for the AGNs and 
star-fo rming central regions selected from the work of I Ho et al.l 
l|l997l) . Solid and dashed curves mark the bounda r y bet ween 
AGNs and Hl l reRio ns defined by iKaufFmann et al. I J2OO3I ') and 
iKewlev et al. I ll200l] 'l. respectively. The filled circles show a large 
sampl e of emission-line SDSS galaxies studied by iThuan et al.l 

il2ninl'i . 


ABUNDANCE AT THE CENTRE OF A 
GALAXY vs. CENTRAL INTERSECT 
ABUNDANCE 


3.1 The data 

It is a widely accepted practice to specify the abun¬ 
dance at the center of a galaxy by the central intersect 
abun dance obtained from the radia l abundance gradient 
(e.g., Vila-Costas fc Edmunds 19921 : Zante^ et al. 1994 
\gnZee et alJ 19981 : Pilvuein et alJ 2004I . 20071 : Gusev et alJ 
201 I among others). The radial abundance gradients in the 
dis ks of nearby la t e-type gala xies were recently determined 
bv IPilvugin et all (l20l4[^15l ). The extrapolation of the ra¬ 
dial abundance gradient to the zero galactocentric distance 
gives the central in tersect abundan ce (0/H)o in galaxies. 
On the other hand, IHo et all (Il997l i obtained emission-line 
spectra of the central regions of many nearby galaxies. This 
provides a possibility to estimate the oxygen abundances 
at the centers of those galaxies. The comparison between 
central oxygen abundance estimated from the spectrum of 
the center of the galaxy with that obtained from the radial 
abundance gradient can tell us something about the chemi¬ 
cal evolution of the central parts of galaxies. 

We compare the central 0/H abundances estimated 
from the spectra oflHo et al.l (Il997li and the central intersect 
abundances obtained by Pilvuein et al.l ll20l4l20f5l 'l. We se¬ 
lected from the sample of Ho et al. (|l997^ mly the galaxies 
with central spectra classified bv lHo et al.l (Il997li as H ll-like 
regions or S eyferts and which are also in the Pilvugin et al.l 
(l20l4 l2015li list. Our selected sample contains 45 objects 


(12 AGNs and 33 star-forming regions) whose emission-line 
intensities and derived oxygen abundances (see below) are 
listed in Table [S] Fig. [4] shows the positions of the 45 sam¬ 
ple o bjects in a BPT diagnostic diagram ([Baldwin et al.l 
H), together wit h a large sample of emission-line SDSS 
galaxies studied by iThuan et al.l (l2010l '). which are plotted 
with cyan (grey) symbols. In this figure we also plotte d the 
boundary line s obta ined by iKauffmann et ai~l (l2003l ) and 
iKewlev et al. I (l200ll ~) that separate the H ii-like objects and 
AGNs in the diagram. 

IPilvugin et al.l (I20l4 l2015l l derived oxygen abundances 
from the published emission-lin e intensities of disk H ll re¬ 
gions trough the C method fsee lPilvugin et ^120121 '). This 
is an empirical method based on the comparison between 
strong emission-line intensities in the spectrum of a target 
HII region and those in a set of reference H ll regions with 
known abundances. The strong emission-lines considered in 
this method are [Oii]A3727, H/3, [O iii]A 5007, [N ii 1A6584 
and [Sii]A6716, A6731. The optical data of lHo et al.l (Il997l ) 
consist of long-slit spectroscopy of the nuclear region (r < 
200 pc) of a large sample of nearby galaxies covering the 
4200-5200A and 6200-6900A spectral ranges with spectral 
resolutions of about 4 and 2.5 A pixel”^, respectively. We 
use again the first relation of SB98 to estimate the central 
abundances of the NLRs of the active galaxies [(0/H)sb98,i]- 
In the case of gala xies with central Hll - like re gions we use 
the Cns method dPilvugin et al.l [ 2 OI 2 I . l2013l i to estimate 
the oxygen abundances [(0/H)cjvsl- The Cns method is 
a variant of the C method applicable when the [Oll]A3727 
emission-lines are not available. 
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Table 3. De-reddened nuclear emission-line fluxes (normalised to the flux of H/3=1.00) taken from lHo et alJ i|l997h and oxygen abundances 
derived by us. The oxygen abundances 12+log(0/H) for the NLRs of AGNs and central star-forming regions are estimated through the 
SB98 calibration (Eq. [TJ and through the Cns method, respectively. The values of 12-|-log(0/H)o are the central intersect abundances 
in the galaxies obtained from the radial abundance gradients. 


Object 

[O iii]A5007 

[N ii]A6584 

[Sii]A6717+A6731 

12-blog(0/H) 

12+log(0/H)o 

NLRs of AGNs 

NGC 1058 

3.64 

3.50 

1.95 

8.58 

8.62 

NGC 1068 

12.10 

2.17 

0.68 

8.54 

8.64 

NGC 2336 

2.84 

5.12 

2.99 

8.74 

8.80 

NGC 3031 

4.07 

6.38 

3.89 

8.78 

8.58 

NGC 3227 

5.92 

3.82 

1.95 

8.62 

8.64 

NGC 3486 

4.50 

3.01 

2.64 

8.56 

8.60 

NGC 4258 

10.20 

2.29 

2.61 

8.51 

8.54 

NGC 4395 

6.23 

0.93 

2.04 

8.34 

8.19 

NGC 4501 

5.15 

6.02 

2.66 

8.80 

8.92 

NGC 4725 

6.68 

3.27 

1.82 

8.59 

8.83 

NGC 5033 

4.59 

6.76 

2.99 

8.85 

8.64 

NGC 5194 

8.09 

8.26 

2.28 

9.01 

8.88 

Nuclear star-forming regions 

NGC 598 

1.57 

0.48 

0.61 

8.35 

8.48 

NGC 783 

0.19 

1.03 

0.55 

8.65 

8.68 

NGC 925 

0.81 

0.63 

0.96 

8.41 

8.48 

NGC 1156 

4.51 

0.28 

0.46 

8.24 

8.16 

NGC 2403 

1.91 

0.79 

1.03 

8.36 

8.48 

NGC 2537 

1.84 

0.43 

0.48 

8.35 

8.35 

NGC 2903 

0.08 

0.97 

0.52 

8.71 

8.82 

NGC 2997 

0.31 

1.00 

0.77 

8.56 

8.80 

NGC 3184 

0.12 

0.94 

0.55 

8.67 

8.66 

NGC 3198 

0.23 

1.19 

0.84 

8.59 

8.60 

NGC 3310 

0.91 

1.88 

0.71 

8.55 

8.37 

NGC 3319 

0.99 

0.46 

1.71 

8.26 

8.50 

NGC 3344 

0.83 

1.20 

1.15 

8.49 

8.72 

NGC 3351 

0.25 

1.31 

0.66 

8.62 

8.82 

NGC 3359 

0.54 

0.88 

1.55 

8.43 

8.40 

NGC 3631 

0.29 

1.23 

0.68 

8.60 

8.71 

NGC 3738 

2.99 

0.26 

0.81 

8.13 

8.10 

NGC 3893 

0.21 

1.05 

0.61 

8.62 

8.73 

NGC 4088 

0.19 

0.91 

0.48 

8.65 

8.71 

NGC 4214 

3.66 

0.19 

0.33 

8.20 

8.20 

NGC 4254 

0.87 

1.37 

0.59 

8.54 

8.77 

NGC 4303 

1.32 

2.23 

1.12 

8.53 

8.78 

NGC 4449 

2.36 

0.40 

0.65 

8.26 

8.26 

NGC 4490 

2.34 

0.71 

1.89 

8.22 

8.29 

NGC 4535 

0.13 

1.17 

0.61 

8.66 

8.71 

NGC 4559 

0.32 

1.20 

1.10 

8.53 

8.53 

NGC 4631 

1.51 

0.69 

0.63 

8.42 

8.39 

NGC 4654 

0.14 

0.77 

0.60 

8.63 

8.66 

NGC 4656 

3.96 

0.14 

0.50 

8.04 

8.06 

NGC 5248 

0.26 

1.48 

0.64 

8.61 

8.64 

NGC 5457 

0.23 

1.08 

0.70 

8.61 

8.71 

NGC 5474 

1.75 

0.37 

0.72 

8.26 

8.19 

NGC 6946 

0.37 

1.81 

0.81 

8.55 

8.72 


3.2 Abundance at the centre vs. central intersect 
abundance 

In the bottom panel of Figure[5]we show the comparison be¬ 
tween the central intersect abundances in the selected galax¬ 
ies (0/H)o obtained from the radial abundance gradients 
and central abundances (0/H)sb98,i and (0/H)c;vg deter- 
m ined from t he sp ectral measurements of the central regions 
of lHo et"^ lll997l i. This figure shows that the (0/H)sb98,i 


abundances in the NLRs are close to or slightly lower than 
the central intersect abundances in the host galaxies. The 
central abundances in galaxies with central H ll-like regions 
show a similar behaviour. Thus the metallicity in the NLRs 
obtained through the relation of SB98 is close to the cen¬ 
tral metallicity of the host galaxy estimated from the ra¬ 
dial abundance gradient. The mean difference (D) between 
the direct and intersect central oxygen abundances (see up- 
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Figure 5. Comparison between central intersect oxygen abun¬ 
dances derived from the radial abundance gradients (0/H)o and 
the central abu ndances determi ned from the spectra of the cen¬ 
tral regions from lHo et alJ l|l997tl through the SB98 calibration for 
the AGNs (filled squares) and through the Cns method for the 
Hii-like regions (open circles). The diagonal line in the bottom 
panel is the one-to-one relation. The upper panel shows the differ¬ 
ences between the spectral and central intersect abundances. The 
hatched area indicates a band of 0.15 dex adoped for the oxygen 
abundance uncertainty (see Section Oi- Solid line shows the lin¬ 
ear regression [D = (—0.25±0.06) X 12-|-log(0/H)o-l-2.10(±0.59)] 
considering all data. 


per panel of Figure [SJ for Hll-like objects is —0.11 ± 0.09 
at high metallicities (12-|-log(0/H)>8.5) and 0.00 ± 0.08 at 
low metallicities (12+log(0/H)<8.5). The mean difference 
for AGNs is —0.01 ± 0.13 at high metallicities. The differ¬ 
ences for the high metallicity AGNs seem to follow a linear 
regression with a slope of —0.32(±0.33) dex“^. The differ¬ 
ences for H ll-like objects follow a regression with a slope of 
—0.45(±0.23) dex“^ at high metallicities, and with a slope 
of —0.23(±0.15) dex“^ at low metallicities. We also per¬ 
formed a linear regression considering all points in Fig. [S] 
due to the small sample we are using. It yields a slope of 
-0.25(±0.06) dex-^ 

The differences are within the uncertainty in the abun¬ 
dance estimations (i.e., compatible with zero). However, the 
existence of a trend in the differences suggests that the di¬ 
rect central abnndance in some high metallicity galaxies can 
be lower th an the central i nterse ct abundance. It should be 
noted that ISanchez et al.l (l2014l i have also found observa¬ 
tional evidence of lower central oxygen abnndances than 
that predicted by the central abundance extrapolation of the 
gradients in a number of spiral galaxies. A possible expla¬ 
nation of the trend of the differences with metallicity could 
be the accretion of metal-poor gas onto the centers of galax¬ 
ies. In the case of the low-metallicity galaxies, where the 
metallicity of the accretion material is similar (or not so 
different) to that of the gas in the central region, the in¬ 
falling gas should not significantly change the local metal¬ 


licity. In contrast, in the case of high-metallicity galaxies this 
low-metallicity infalling gas would dilute the heavy element 
content of the gas in the central region. Therefore this effect 
would be more relevant in high-metallicity galaxies than in 
low-metallicity ones. 

Observational evidence of the presence of gas flows 
from the outer parts (with low-metal content material) to 
the centre of the galactic disc (high metal content) has 
been found for isolated barred galaxies (|Marti^_&_^^ 1994; 


Zaritskv et al. Ill994l) and i nteracting ones ( [Rosa et al 


20141 : 


Ellison et al. 201 ll . [20101 : iKewlev et al.l l2010l i. Likewise 
there are several kinematical studies based on optical and 
infrared integral held spectroscopy that support the scenario 
where gas is infahing toward the c e ntral region of AGNs (e.g . 
Schnorr-Muller et al.l 12(314 l201ll: i Riffel et lj] l2013bl . 1200^ : 


Storchi-Bergmann et al.l 2007 : Fathi et al. 2006I L Moreover, 
Rupke et al.l 1 2010l 'l ^rformed numerical simulations of 
galaxy mergers and found that the central underabundances 
observed in this kind of systems could be accounted for by 
a radial infall of low-metallicity material coming from the 
outskirts of both galaxies involved. 


3.3 Is there an extraordinary chemical 
enrichment of the NLRs? 

It was found by IPilvugin et al.l (l2006l . l2007l i that there is 
an upper limit to the oxygen abundances in galaxies, i.e., 
there is a maximum attainable oxygen abundance. These au¬ 
thors found that this maximum value of oxygen abundance 
in galaxies is about twice the solar abundance (adopting 
the solar oxygen abundan ce to be 12 + log(0/H)Q = 8.69; 
lAllende Prieto et al.ll200i] '). 

The a bundance s in t he NLRs of AGNs taken from 
the list of et al.l (1 19971) are in the range 8.5 < 12-1- 
log(0/H) < 9.0 or 0.6 < Z/Zq < 2, which does 
not exceed the maximum attainable abundance for galax¬ 
ies. This suggests that there is no an extraordinary chem¬ 
ical enrichment o f the NLRs. Similar AGN abundances 
were obtained by iMatsuoka et al.l (l2009l ) from the analy¬ 
sis of the CivA1549/HeiiA1640 - Giii]A1909/CivA1549 di¬ 
agram calibrated in terms of gas abundance through pho¬ 
toion iz ation model s using the Cloudy code (iFerland et al.l 
Il996l) . iDors etld] (l2014h also obtained similar abun¬ 
dances using a new index (C 43 = log[(C lvA1549 -I- 
C III] A1909)/HeIIA1640]) defined by them as a metallicity 
indicator for AGNs. 

However, it should be noted that we found metallic¬ 
ities larger than twice solar up to ZjZ^ ~ 6.5 [12 -|- 
log(0/H)sB98,i ~ 9.5 (see Fig. [T]) for five AGNs from Table 
[T] Similarly extra high metallicities for some NLRs were ob¬ 
tained by SB98. Further study of the galaxies with possible 
extra high metallicity AGNs should be carried out (in par¬ 
ticular, radial abundance distributions and central intersect 
abundances should be obtained in those galaxies) in order 
to be able to draw solid conclusions about the upper limit 
to the oxygen abundances in NLRs. 


4 SUMMARY AND CONCLUSIONS 

We compiled from the literature a sample of spectra of 
narrow-line regions of active galactic nuclei with avail- 
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able optical emission lines: [Oll]A3727, [Olll]A4363, H/3, 
[Oiii]A5007, Hq, [Nii]A6584, [Sii]A 6717, and [Sii]A6731. We 
estimated the oxygen abundances in those NLRs through the 
classic Te me thod and through the strong-l ine method (the 
calibration of IStorchi-Bergmann et aDll998l l . We found that 
the abundances determined through the Te method are lower 
by up to ~2 dex than the abundances estimated through the 
calibration of SB98, i.e., we confirmed the existence of the 
so-called “temperature problem” in AGNs. 

We also considered a second sample of galaxies for which 
the emiss ion-line s p ectra of the central regions were mea¬ 
sured by IHo et all lll997tl and the c e ntral inters ect abun¬ 
dances were found by Pilvugin et alJ (l2014l . [201^ . The di¬ 
rect central abundances of the AGNs and Star-forming re¬ 
gions i n those galaxies were estimate d through the calibra¬ 
tion of IStorchi-Bergmann et ^ (Il998ll and through the Cjvs 
method I Pilvugin et al. I l20lll2013lh respectively. We found 
that the abundances of the NLRs and Hll-like objects esti¬ 
mated from the direct spectral measurements are close to or 
slightly lower than the central intersect abundances obtained 
from the radial abundance gradients. This may suggest that 
the infall of the low-metallicity gas onto the centers of the 
galaxies can take place in some galaxies where the central 
abundance estimated from the direct spectral measurements 
is lower than the central intersect abundance. 

The abundances in the NLRs of the AGNs in our sam¬ 
ples do not suggest that there is an extraordinary chemical 
enrichment of the narrow-line regions. There are only a few 
AGNs with oxygen abundances higher than the maximum 
attainable abundance for g alaxies (~2 times the solar value; 
IPilvugin et al.|[2006l . 12007 1. Additional investigations of the 
galaxies with possible extra high metallicity AGNs are nec¬ 
essary to draw reliable conclusions on the upper limits of 
the oxygen abundances in NLRs. 
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